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We analyze distinct sources of spin-dependent energy level shifts and their impact on the tunneling 
magnetoresistance (TMR) of interacting quantum dots coupled to collinearly polarized ferromagnetic 
leads. Level shifts due to virtual charge fluctuations can be quantitatively evaluated within a 
diagrammatic representation of our transport theory. The theory is valid for multilevel quantum 
dot systems and we exemplarily apply it to carbon nanotube quantum dots, where we show that 
the presence of many levels can qualitatively influence the TMR effect. 
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Recent transport experiments on quantum dots cou- 
pled to ferromagnetic leads have demonstrated the ex- 
istence of spin-dependent energy level shifts, through 
the observation of negative tunnel magneto-resistance 
(TMR) effects in the single electron tunneling regime [1] , 
and spin splitting in the Kondo regime [2-5] . So far quali- 
tative different approaches for explaining the origin of the 
underlying shifts co-exist. For example, negative TMR 
data from CNT measurements have been fitted with a 
model relying on spin dependent interfacial phase shifts 
[6] picked up by the wave function during multiple reflec- 
tions at a spin-active interface [7] . The concept is related 
to that of spin-mixing conductance [8] and, because it 
only depends on the properties of the spin-active barrier 
region, is only weakly gate dependent and present both 
in interacting and in non-interacting systems [9] . In con- 
trast, experiments on CNTs [3, 5] and InAs nanowires [4] 
in the Kondo regime have demonstrated a combination 
of gate dependent and gate independent contributions to 
the energy level shifts. Both can be explained in terms 
of charge fluctuations [10], whereby electron-electron in- 
teractions are responsible for the logarithmic gate depen- 
dence [11], while a Stoner splitting of the energy bands of 
the magnetically polarized leads accounts for the almost 
gate independent part [12]. 

While the effects of the energy level shifts in the Kondo 
regime are by now well understood, a thorough under- 
standing of their influence on the TMR phenomenon in 
interacting quantum dots is still missing. The negative 
TMR data [1] have been satisfactorily fitted in terms 
of a generalized Anderson model already including gate- 
independent level shifts. In [13] a reflection Hamiltonian 
was included to account for reflection processes at the 
interface in a second order sequential tunneling theory. 

In this work, we specifically address the TMR phe- 
nomenon and discuss different level-shift induced mech- 
anisms yielding negative TMR. In particular we discuss 
two types of spin-dependent level shifts that arise intrin- 




FIG. 1: (a) Transport set-up of a carbon nanotube (CNT) 
quantum dot with ferromagnetic leads. If an electron enters 
the barrier at the tube ends, it can either tunnel to the lead 
or being reflected at the interface. Tunneling processes in the 
contact region are described by even powers of the tunneling 
Hamiltonian Ht- (b)/(c) Conductances Gp,Ga versus gate 
voltage for two situations exhibiting negative TMR. (d) The 
negative TMR mechanism related to (c), where thick/thin 
lines describe processes that are favoured/disfavoured being 
associated to majority/minority spins. 

sically from tunneling induced renormalization: strongly 
gate-dependent ones from a net difference of majority and 
minority density of states at the Fermi level, and largely 
gate- independent ones from a Stoner shift of the majority 
and minority bands in the leads. We show how to include 
these effects within a diagrammatic approach to the re- 
duced density matrix of the nanosystem [14] for generic 
multilevel quantum dots, and reproduce the results of 
[12] in the case of a simple Anderson model. Finally, 
we analyze the TMR of a CNT quantum dot, and show 
that, due to the multilevel spectrum, intrinsic contribu- 
tions only can yield a marked gate voltage dependence 
and a TMR which can indeed become negative. 

Mechanisms of negative TMR. The basic mechanism 
underlying a negative TMR is the presence of an effective 
generalized Zeeman field h p / A (V s ) accounting for both 
extrinsic (stray fields, applied magnetic field) and intrin- 
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sic sources of spin-dependent level shifts. Here (P/A) 
refers to contacts with parallel/antiparall magnetization, 
see Fig. 1. Intuitively, one would expect that the con- 
ductance for contacts polarized in parallel (Gp) is larger 
than the one for the anti-parallel case {Ga), yielding- 
positive values of TMR := (G P /G A ) - 1. Neverthe- 
less, there are at least two different mechanisms which 
can lead to a negative TMR, as sketched in Figs. 1(b), 
(c), respectively. Here we consider quantum dots in the 
single electron tunneling regime and look at the linear 
conductance as a function of gate voltage V g . For suffi- 
ciently small couplings the temperature determines the 
width of the conductance peaks, while the peak position 
signals at which value of V g the two many-body ground 
state configurations with N and N + 1 electrons involved 
in the transition have the same electrochemical poten- 
tial: /j,(N,V g ) = p(N + 1, V g ). Because the effective 
Zeeman field leads to corrections to the energy differ- 
ence Ei, — E a := E oa associated to the transition from a 
many-body state a to a many-body state b, this in turn 
modifies the position of the conductance peaks. 

First, and most obviously, there can be negative TMR 
if there is a noticeable offset in the conductance peak 
positions in parallel and anti-parallel configuration, see 
Fig. 1(b). This requires \h p — h A \ of the order of the 
width kpT of the conductance peaks. Secondly, however, 
negative TMR can also arise for \h p \ = \h A \ « kpT, see 
Fig. 1(c). This is because the effective magnetic field 
responsible for the effective Zeeman shift also removes 
the spin-degeneracy of the ground states by favoring the 
states with maximum total spin. This situation is illus- 
trated in Fig. 1(d) for the case of a CNT quantum dot, 
for which we consider the CNT Hamiltonian Eq. (2). The 
four-electron shells of the dot are filled sequentially upon 
sweeping the gate voltage and only a single shell needs 
to be considered at a time. The case of orbital (Ae = 0) 
and spin degeneracy is illustrated for the n-th shell in 
Table I. For the An o An + 1 and An + 1 o An + 2 tran- 
sitions and h p l A ^ only spin-up electrons are required. 
Hence, due to the larger density of spin up electrons, 
peak heights for parallel contact polarization will exceed 
those for the anti-parallel case. On the other hand, for 
the An + 2 o An + 3 and An + 3 <-> An + A transitions spin- 
down electrons have to be transferred through the dot. In 
this latter case a configuration with anti-parallel contact 
polarization might be favored leading to Ga > Gp. 

The model Hamiltonian For a quantitative descrip- 
tion we consider the transport setup of Fig. 1. In the 
limit of weak coupling to the leads it can be described by 
the total Hamiltonian 

Ki A = H + H P J? + Hi + H T , (1) 

l=s,d 

where H = H^ — eaVgN comprises the Hamiltonian Hq of 
the isolated quantum dot and the effects of a gate voltage 
(a is a conversion factor of the order of one). In the 



case of an armchair nanotube quantum dot of medium- 
to-large radius far from half-filling it reads [15] 

H = l -E c N 2 + X - Y,(£oN fa + rAe)N fa + Hp, (2) 

f ,CT 

where e = Jivf/ttL, with v-p the Fermi velocity and 
L the CNT length, is the level spacing; Ae is the or- 
bital mismatch, E c is the charging energy and Hp ac- 
counts for bosonic charge and spin excitations. The num- 
ber of electrons in orbital band f with spin a is N fa 
and the total number is N = ^2 f<T Nf a . The ground- 
states | N + f, N +l , 7V_ t , 7V_ t ) of shell n have An, An + 1, 
An + 2 and An + 3 electrons and can be characterized 
in terms of the excess spins in band {f, a} with respect 
to the case of equally filled bands: | n,n,n,n) :=| •,•), 
n + 1, n,n, n) :=|t, ■), or | n + l,n + 1, n,n) :=|t4-, •)> 
etc., see Table 1. For medium-to- large tube radius and 
far from the charge neutrality point spin-orbit coupling 
and exchange effects are not relevant. They can become 
of interest for large curvatures and will be investigated 
elsewhere. H^J^ accounts for external, gate independent 
sources of level splitting. 

The leads are described by H t = Y, aq ei q cj aq ci aq , with 
Q aq annihilating an electron of energy e; 8 and of spin a 
in lead I. We assume a constant density of states in lead 
I, D la , which determines the leads polarizations P; at the 
Fermi level according to Pi = (Di-f — Di^)/{D^ + Di±). 
On the other hand, in order to account for a Stoner 
splitting Ag t , the range of available energies in the 
leads is given by [12] -{W + oA St /2) < w < W - 
crAgt/2, where W is the bandwidth at zero Stoner split- 
ting. Finally the perturbative contribution is H T = 

J2i a fd ? ' r (Ti(r)i>^(r)4> la (r) + h.c.^j , allowing for tunnel- 
ing between CNT and leads. Here T;(r) is the tunnel 
coupling, the CNT bulk electron operator as given 

in Ref. [15], and 4> la (r) = J deD la (e) £,| e 4> iq (r)c laq the 
lead electron operator with 4>i q {r) denoting the corre- 
sponding single-particle wave function. 

Transport theory. The current, as any other observ- 
able of the transport set-up, can be calculated by a trace 
over the associated operator multiplied by the reduced 
density matrix p{t) = Tr; {p tot (t)}. Being obtained from 
the density matrix p tot of the total set-up upon tracing 
out the lead degrees of freedom, p(t) stores the full infor- 
mation about the state of the dot in the presence of the 
tunnel coupling to the leads. The time evolution of p(t) 
follows the Liouville equation, and it can for the station- 
ary state (p(t) — 0) be cast into the form, see e.g. [16], 

= -t X) 5 "b Sa'V (Ea ~ K> ) Paa> + ^ K W Paa' , (3) 
aa' aa' 

taking matrix elements with respect to the many-body 
eigenstates of the CNT: p aa > := (a\p\a'} and E a := 
(a\H\a). Furthermore K$ := (b\[K(a)(a'\]\b') , with the 
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filling 


ground state(s) 


N = An 


It) 


N = An + 1 


ItrJ.ht) 


,|4,-),h4) 


N = An + 2 


|t,t> ,lT4,-},lT,4},|-,t4},l4,t},l4,4) 


N = An + 3 


It4,t),lt,t4) ,IT4,4),I4,T4) 



TABLE I: Ground states of a CNT in shell n accounting for 
spin a = {t,4} an d orbital degrees of freedom r = {+/—}, 
leading to four-electron shells. The number of electrons in 
band f with spin a is Nf a and in the dot N = ^ f Nf a . 
With h p/A / 0, only the highlighted states are ground states. 



kernel superoperator K arising from the perturbation, 
i.e., the tunnel coupling between quantum dot and leads. 
Depending on up to which order 2n in the tunnel cou- 
pling K is calculated, one takes into account effects from 
In correlated tunneling events. The most involved part 
of a perturbative transport calculation is the determina- 
tion of the kernel K. In the supplemental material we 
identify specific terms in all orders of the perturbation 
series as charge fluctuation processes, which sum up to a 
Taylor series yielding the intrinsic level renormalization 
E a b — > Eab + hf b t to the energy difference E a b, where 



-£K( c >)f/ 



du>- 



dto 



f(u)Di a (u) 
- f3E cb - pfleVi 

(1-/H)A.H 



(4) 



We notice that Eq. (4) generalizes to complex quantum 
dots the analysis results for the energy level shifts ob- 
tained for the single impurity Anderson model [10-12]. 
It is convenient to rewrite h\® t = h b a (V g ) + Ah ba {A st ), 
where A/i 6a (As t = 0) = 0. As shown in the supple- 
mental material, the first term depends logarithmically 
on the gate voltage and the interaction strength, while 
Ah ba depends on the finite bandwidth W a = W-a A st /2 
of D[ a (ijj) and will be almost gate and interaction inde- 
pendent and thus contribute in the same way as contri- 
butions h^lf induced by H^J t A . 

Notice that for quantum dots with PdNi contacts, 
given a tunneling coupling T ~ 1 meV, one gets Ah ba 
of the order of few hundreds jieV [5], This is the same 
order of magnitude as that of the effective constant field 
used to fit the TMR data in Ref. [1] with the interface 
phase shift model [6] ( T = J2i T i and r * ~ A I T i | 2 )- A 
full account of the interface phase shifts might be neces- 
sary for large transparencies and goes beyond the scope 
of our tunneling Hamiltonian. However, in the regime 
where the tunnelling model is valid (Coulomb blockade 
and Kondo regime), tunneling induced shifts should suf- 
fice to explain the experiments. 

Linear conductance results. We consider the linear 
bias regime for spin polarized transport across a CNT 



without band-offset, cf. Table 1. Fig. 2 shows the con- 
ductances for parallel (pink solid) and anti-parallel (blue 
dashed) lead magnetizations under the influence of the 
distinct types of level shifts, along with the correspond- 
ing TMR curves. Throughout this work we assume equal 
lead polarizations P s = Pa = P and a coupling asymme- 
try 7 = r ' s/Ya ~ |T s | 2 /|Td| 2 , with the actual values of 7 
given at each plot. 

In the absence of any Zeeman shift, Fig. 2(a) shows the 
typical conductance peak patterns expected for the de- 
generacies according to Tab. I, for both a second order 
(sequential tunneling) and a fourth order (cotunneling, 
pair tunneling, single charge fluctuations etc.) [16] trun- 
cation in the calculation of the kernel K. Second or- 
der theory predicts a constant TMR value [13], but this 
changes upon inclusion of higher order effects. To fourth 
order in perturbation theory the TMR exhibits an os- 
cillatory gate voltage dependence, albeit the variation 
is small [17]. To visualize the TMR mechanism from 
Fig. 1(c) and (d), we include in Fig. 2(b) merely a con- 
stant shift from an external source, \hg Xt \ = \h£ xt \ > kgT 
(this can arise from the external field which controls the 
direction of contact polarization in TMR measurements, 
when for the measurement of Gp the field is swept back 
to the point where Ga has been measured before). As 
explained above, among the formerly degenerate lowest 
lying states, the ones with maximum spin-projection S z 
are selected as ground states, as highlighted in Tab. I. At 
the first two peaks transport is mediated by \ -, at the 
last two peaks by 4 - electrons. For parallel magnetiza- 
tion of the leads this breaks the mirror symmetry of the 
conductance by enhancing the first two peaks while sup- 
pressing the last two. In Fig. 2(c) we still neglect intrisic 
shifts, but assume \h^ xt \ > \h^ xt \ (> kgT, e.g. due to 
stray fields) to bring also the first negative TMR mecha- 
nism, Fig. 1(b), into play. It overlays the effects observed 
in Fig. 2(b). The shift in the peak positions of parallel 
and anti-parallel conductance peaks enforces a change of 
the TMR from positive to negative value from the first 
peak on. Because there is no significant change in the 
curves of Fig. 2(b) and 2(c) upon inclusion of the fourth 
order effects, in the remaining figures we only show sec- 
ond order curves. In Fig. 2(d) we include the finite level 
shifts which arise intrinsically from our transport theory, 
both for Ast = and Ast = 0.2W. In the latter case 
we find that inclusion of excited states (here up to ex- 
citation energies of 2cq) can even yield negative TMR. 
Interestingly, the intrinsic gate dependent contribution 
distinguishes itself from all former contributions by pre- 
serving the mirror symmetry in the conductances and 
the TMR curve, as seen for Ag t = 0. The reason is 
that a given electron-hole symmetry in the tunneling re- 
sults in an electron-hole symmetry for any charge fluc- 
tuation (see supplemental material). This is no longer 
true if Agt 7^ 0, where we observe a behaviour similar 
to that of a constant external shift but of opposite sign 
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FIG. 2: Parallel (Gp) and anti-parallel (Ga) conductance along with the resulting tunnelling magneto-resistance (TMR) for a 
CNT of 500 nm length (eo = 3.35 meV) and charging energy E c = 6.7 meV. The thermal energy was set to fcgT = 0.3 meV, the 
lead polarization to P = 0.4, the tunnel broadening to KT S = 3 /^eV and the bandwidth to W = 3eV. (a) For a full fourth order 
calculation a mirror-symmetric TMR slightly oscillating around a value of 20% is obtained, (b) An equal splitting h^xt — h^ xt = 



0.4meV causes gate asymmetry and negative values of the TMR by the mechanism Fig. 1(c). (c) For h e3:t = 2h B 



0.8 meV, 



the regularity of the curve is broken and the negative TMR mechanism Fig. 1(b) comes into play, (d) TMR for A st = 0, 0.2W 
and conductance for Ast = 0. Including the influence of the excited states (with cut-offset to eo and 2eo) negative TMR can be 
reached. With Ast = electron-hole-symmetry is preserved, for Ast ^ is broken, (e) This is confirmed by direct investigation 
of the intrinsic level shifts for some electron-hole symmetric transitions in parallel configuration, (f) Combining intrinsic shifts 
with external splitting. An asymmetric coupling to the leads can enlarge (7 < 1) or diminish (7 > 1) the TMR effect. 



[5]. These statements are confirmed by Fig. 2(e) which 
shows the shifts for exemplary electron-hole symmetric 
transitions. Ag t 7^ induces a gate constant positive 
(negative) shift, for f (I) - mediated transitions, which 
breaks the electron-hole symmetry. It is nicely visible 
that, due to the orbital degeneracy of the CNT, the shifts 
differ from peak to peak. Finally, Fig. 2(f) combines 
the impact of the intrinsic shifts with an equal extrin- 
sic splitting \he Xt \ = \h^ xt \. The intrinsic effects suffice 
to change the TMR curve Fig. 2(b) to a shape observed 
in experiment [1], though due to our limitation to the 
weak coupling regime, quantitative agreement cannot be 
expected. An asymmetric coupling is found to affect the 
curve quantitatively, while the qualitative shape is re- 
tained. Even for the very small values of T s used here 
to justify lowest order perturbation theory, the difference 
between Figs. 2(b) and (f) reveals a marked influence of 
the intrinsic, tunneling induced, level shifts. 

Summary We analyzed the impact of different kinds 
of effective Zeeman shifts in magnetically coupled mul- 
tilevel quantum dots, obtaining a characteristic gate de- 
pendence and the possibility of negative TMR. In par- 
ticular, we have provided a systematic way of including 
the important effects of tunneling induced level shifts in 



a transport calculation. In general, a TMR signal will 
be influenced by many parameters relevant to the given 
device. Nevertheless, following the lines of the analysis 
given here for a CNT, it should be possible to disentangle 
the importance of the various contributions. 

We acknowledge fruitful discussions with C. Strunk, S. 
Pfaller and J. Hauptmann as well as support of the DFG 
under the program SFB 689. 
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In the main part of our Rapid Communication, we have 
introduced with Eq. (2) the Liouville equation for the re- 
duced density matrix p(t) = TV; {ptot(^)} m the station- 
ary state (p(t) = 0), 

= ~i & ab <W {E a ~ E a > ) p aa > + ^ K$ Paa > ■ 
aa' aa' 

The solution to this equation describes the stationary 
state of the quantum-dot system in the presence of the 
leads. Thereby, the main task is the determination of 
the kernel K. A fully diagrammatic representation of K 
has first been proposed in [1]. For example the second 
order kernel {K( 2 >)%£ is shown in Fig. 1(a). It describes 
a full tunnel event \a) — > b) , represented by the solid line 
connecting upper and lower parts of the diagram contour. 
Recently, it has been recognized [2] that certain fourth or- 
der diagrams can be related to charge fluctuations in the 
initial or final state during a tunnel event. Here we iden- 
tify a specific class of these charge fluctuation diagrams 
which, summed up in all orders, yield a whole Taylor se- 
ries and therewith an actual level shift. Fig. 1(b) shows 
a diagram of order 2n + 2, contributing to (jf ( 2 ™+ 2 ))a<^ 
which dresses Fig. 1(a) by k charge fluctuations in the fi- 
nal state b (k "bubbles" on the upper part of the contour) 
and n — k charge fluctuations in the initial state a (n — k 
"bubbles" on the lower part of the contour). We thereby 



look at fluctuations isolated, in the sense of separated in 
time, from each other: each bubble must start and end 
at consecutive times tj, r^+i (1 < i,j < 2n — 1). The 




FIG. 1: Structure of: (a) a sequential tunneling diagram, 
(b) a level renormalization diagram for a sequential tunneling 
event. Each "bubble" describes a charge fluctuation. We 
consider k charge fluctuations (to states Cj) in the final state 
b and n — k charge fluctuations (to states Cj) in the initial 
state a. We have indicated the time order, as well as the 
intermediate states on the contours. 

electron transfer event is initialized at the earliest time 
r = To, and ending at the latest time T2 n +i = t. This 
gives (?) possibilities for the time ordering of the bubbles 
among the upper and lower part of the contour. Sum- 
ming all those plus their hermitian conjugates, the total 
contribution is, as obtained by applying diagrammatic 
rules [1], 
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Here Cj, c! i serve as placeholders for the possible inter- a are set by fixing the initial and final states a and b. 
mediate states; := / p (/3o; + /3eVj), p = ±, where Likewise the values of the spin indices aj, o~i and of the 



we defined f + {u)) = /(w) = 1/(1 + e"), denoting the 
Fermi function, while f~(to) = f{—uj). The lead indices 
l,h,lj € {s,d} are summed over. The values of p and 



Fermi function labels pj, pi depend on the intermediate 
states Cj, d i that can be reached from b, a, respectively, 
and are summed over. Here, E^a stands for the energy 
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difference E b — E a between the many body states of the 
quantum dot and the lead energies are integrated over. 
Finally, there are tunnel matrix elements characterizing 
transitions on the dot, T^ (+) (&,a) ~ T;(r;)(£#S +) (n)|a) 
(?i characterizes the position of lead I). A multitude of 
terms emerges from the expression Eq. (1) when apply- 
ing the decomposition 



lim / 



Df p H 
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(J' denotes a principal part integration) for all fractions 
and expanding the product. We want to focus on a cer- 
tain one, namely the combination where the fraction con- 
taining merely cu, 
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has been replaced by the delta function, and all others 

by their principal parts. 

In terms of f^(uj) = a^b/H [t reads: 
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Exhausting all possibilities of choosing the states Cj, c[ 
as well as setting k — 1 • • • n, it is clear that one generates 
all kinds of terms appearing in an expansion of the power 
n of a quantity h\® t defined as 



^ = £(£|Tf CT M)| 2 /' d -- 



■ pE cb - ppeVi 



2 *f(l3(E ba + 



where c and c' run over all states connected to |6), respec 
tively | a) via a charge fluctuation. In total, we obtain the 
nth term ^^fM^E^+peVt)) \Tf a (b,a)\ 
in the Taylor expansion of a Fermi function n 
eVi + h^)) \Tf a {b, a)\ 2 D la . So, effectively, the considered 
contribution of the initial and final state charge fluctua- 
tions renormalizes any energy difference E ba to E ba +h^ t . 
From comparison with Eq. (1), we can extract a renor- 
malization to a manybody energy E b by looking at all 
the possible fluctuations on the upper contour: 



puj 



E, 



cb 



(3) 



where the values of p and a are fixed by the diagrammatic 
rules [1] once the state c is assigned. For the case of the 
single impurity Anderson model this expression reduces 
to the result for the level shift obtained in perturbation 
theory [3, 4]. Hence Eqs. (2), (3) provide a diagrammatic 
interpretation of these results. For a multilevel quantum 
dot, the summation over the intermediate virtual states c 
implies that the amplitude of the level shifts is expected 
to vary from a resonance peak to the other, as also shown 
in Fig. 2(e) of the main part of the manuscript. Notice 
also that for the CNT used in our calculations, the sum- 
mation over the virtual states also include states with 
bosonic excitations as described by the bosonic Hamilto- 
nian Hb in Eq. (2) of the main part of manuscript. 

Having identified the underlying diagrams, it is now 
clear how to include the effects of tunneling induced level 
shifts quite generally in a diagrammatic evaluation of the 
kernel K. As temperature goes to zero, these diagrams 
diverge logarithmically at the charge-degeneracy points, 
and hence give a parametric selection of these specific di- 
agrams. In general, diagrams with multiple single charge 
excitations are expected to give the dominant contribu- 
tion to the charge fluctuations in the presence of interac- 
tions. This infinite class of diagrams is seen to correspond 
to a perturbative renormalization of the dot spectrum, an 
effect which is indeed expected on physical grounds. In 
conventional Feynman diagrammatics, these terms arise 
from summing up the 2nd order dot electron self-energy 
in the Dyson equation. 

In order to better analyze the contribution of majority 
and minority spins we consider the case of a flat band 
including a Stoner splitting Ag t . The density of states 
assumes the form [5, 6] 

D la {u) = D lr7 8{oj+W+aA St /2)d{W-aA S t/2-uj), (4) 

where Di a = Di(l + a Pi) is the density of states at the 
Fermi level, and P; the polarization of contact I. 
Notice that in order to retain holomorphic functions in 
Eq. (2), one can make use of the Fermi function repre- 
sentation of the step functions, 6(±x) — lim^oo f^(j3x). 
Using then the residue calculus [7] , one obtains 



where the upper/lower boundary ±W — crAst/2 of the 
integration enters via 



h 1nt = £ X>„|7£(c,&)| ^_((lE cb +pPeV) 
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where $^0) is the digamma function. An approximation 
of this result is obtained elegantly when using Eq. (4) to 
split the integration range of the integrals in Eq. (2) as 



,W-A st /2 ^.-lV+Ast/2 r 
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dcu = duj + dcu , 
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and noting that in the region at the top/bottom of the 
band is f(uj) = 0/1. Here we get the result 
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where IR^ = W — a Agt/2 is the finite bandwidth for 
the integration over the lead energies. The result is an 
approximation of Eq. (5), as can be seen when applying 
ln(A/_B) « Re*(0.5 + iA) - Re*(0.5 + iB) in Eq. (6). 

Obviously, the last two terms in Eqs. (5), (6) vanish 
if Ast = 0. Moreover, for E c b, E ac > <C W a they are very 
weakly dependent on the gate voltage and polarization, 
but dependent on the saturation magnetization via the 
Stoner splitting. Following [8], with tunneling coupling 
T ~ lmeV [9], this contribution can yield for PdNi con- 
tacts the gate-independent part of the level shift to be 
in the order of few hundreds of /ieV. This is in the same 
order as the constant shift used to fit the Kondo data in 
[8, 10] but also the TMR data [11] in [12]. In contrast, the 
first two terms depend strongly on gate voltage, interac- 
tion and polarization. They correspond to the contribu- 
tion to the level shift first proposed in [4] . This allow us 
to recast the expression for in the form 

h% t = h b a + Ah ba , (7) 

where Ah ab is the contribution due to Ast- 



If the states a, b, c, c' have electron-hole symmetry 
partners a*, c*, c^, there exists the mirror transition 
b± — > a*, with a* now being the final, the initial state. 
Thus the contributions from the gate dependent part of 
the charge fluctuations to is negative, to positive, 
such that h b * a * = —ftp . This preserves the mirror sym- 
metry of all curves when Ast = 0. 

Importantly, this concept generalizes to higher orders: 
virtual charge fluctuations during inelastic cotunneling 
(appearing first in sixth order) give a renormalization 
of the inelastic cotunneling threshold, as already exper- 
imentally observed in [13]. In the present work the con- 
stant part of the level shift has been implemented in the 
numerical calculation of the 4th order contribution to the 
TMR shown in Figs. 2(b) and 2(c) of the main part of 
this Rapid Communication. However, in the considered 
weak coupling regime there is no significant change in 
the curves of Fig. 2(b) and 2(c) upon inclusion of the 
fourth order effects. Because the numerical implemen- 
tation of the gate dependent contribution h b a implies a 
summation over a large number of states for a CNT quan- 
tum dot and no significative changes are expected, in the 
remaining figures Fig. 2(d) - (f ) only the results of a sec- 
ond order calculation with all level shifts fully included 
are shown. 
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